On the Non-equilibrium Phase Transition in Evaporation-Deposition Models 
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We study a system of diffusing-aggregating particles with deposition and evaporation of monomers. 
By combining theoretical and numerical methods, we establish a clearer understanding of the non- 
equilibrium phase transition known to occur in such systems. The transition is between a growing 
phase in which the total mass increases for all time and a non-growing phase in which the total 
mass is bounded. In addition to deriving rigorous bounds on the position of the transition point, we 
show that the growing phase is in the same universality class as diffusion-aggregation models with 
deposition but no evaporation. In this regime, the flux of mass in mass space becomes asymptot- 
ically constant (as a function of mass) at large times. The magnitude of this flux depends on the 
evaporation rate but the fact that it is asymptotically constant does not. The associated constant 
flux relation exactly determines the scaling of the two-point mass correlation function with mass 
in all dimensions while higher-order mass correlation functions exhibit nonlinear multi-scaling in 
dimension less than 2. If the deposition rate is below some critical value, a different stationary state 
is reached at large times characterised by a global balance between evaporation and deposition with 
a scale-by-scale balance between the mass fluxes due to aggregation and evaporation. Both the mass 
distribution and the flux decay exponentially in this regime. Finally, we develop a scaling theory 
of the model near the critical point, which yields non-trivial scaling laws for the critical two-point 
mass correlation function with mass. These results are well supported by numerical measurements. 



I. INTRODUCTION 



Simple lattice-based particle systems undergoing dif- 
fusion and irreversible aggregation provide a useful ap- 
proach to modelling various phenomena in physics such 
as polymer growth kinetics, atmospheric aerosol forma- 
tion and thin film growth on substrates. For a review, 
see [l[ and the references therein. Through various map- 
pings and analogies, diffusion-aggregation models can be 
connected to a quite diverse range of topics in statistical 
physics such as the geometry of river networks Q, self- 
organised criticality @ and force fluctuations in granu- 
lar bead packs It has also been long recognised [|[ 
that such models, in the presence of an external source 
of monomers, also provide excellent examples of systems 
which reach a non-equilibrium stationary state with a 
flux of a conserved quantity (in this case, mass) and ex- 
hibit scaling behaviour far from equilibrium. The sim- 
plest such model, and the basis for much of what follows, 
is Scheidegger's river network model @ or the Takayasu 
model 0, H[ in which particles diffuse with mass indepen- 
dent diffusion rate, coagulate on contact conserving mass, 
with input of particles from outside. Much is known 
about statistics of this model (see [1, [Io| for some more 
recent developments). In particular, it always reaches a 
stationary state and has critical dimension, d c = 2. In 
dimension d > 2 the Smoluchowski mean field theory is 
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applicable and the stationary cluster mass distribution 
scales as to -3 / 2 whereas for d < 2 diffusive fluctuations 
become essential and the scaling of the cluster mass dis- 

2d + 2 

tribution is modified to m d + 2 . In all cases the sta- 
tionary state is characterised by a constant flux of mass 
from small to large masses. As a consequence, the scal- 
ing exponent of the 2-point mass correlation function re- 
sponsible for the mass flux (to be defined later) can be 
determined exactly in all dimension, a result we refer to 
as a Constant Flux Relation (CFR) [III- 

In this paper we are interested in studying what hap- 
pens when monomers are allowed to evaporate from clus- 
ters. By allowing evaporation, mass is no longer strictly 
conserved. Several questions then arise which we attempt 
to answer in this paper. Does evaporation change the 
character of the stationary state? In particular, what 
happens to the CFR? The effect of evaporation was con- 
sidered in fl2| in a closed system without injection. In 
that study, evaporating monomers remained in the sys- 
tem acting as an effective source so that the total mass 
was exactly conserved. The effect of evaporation in an 
open system (in the sense that evaporating monomers 
leave the system) with an external source of monomers 
was studied in [l3|. As a result of these studies, it was 
understood that the system exhibits a non-equilibrium 
phase transition as the rate of evaporation of monomers 
is varied. For weak evaporation, the constant flux sta- 
tionary state persists, albeit with a modified mass flux. 
In hj we gathered numerical evidence that the CFR 
scaling exponent was insensitive to evaporation in this 
regime. In what follows, we make this observation pre- 
cise. This regime does not, however, persist for all values 
of the evaporation rate. Above a certain critical value 
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of the evaporation rate, it becomes impossible to gener- 
ate large clusters and the mass flux decays exponentially 
with cluster mass. In this regime, we show that the CFR 
is replaced by a scale-by-scale balance between the ag- 
gregation and evaporation fluxes. In this paper we refer 
to these two regimes as the growing phase and the expo- 
nential phase respectively. 

The connection between the flux of a conserved quan- 
tity in a non-equilibrium steady state and universal as- 
pects of the underlying statistics has been an important 
theme in the theory of fluid turbulence (see [l5l - [T7j for 
a review). Kolmogorov's famous 1941 paper identified 
the flux of energy through scales of a turbulent flow as 
the single most important feature of a turbulent steady 
state. In particular, for stationary, isotropic turbulence, 
the constant flux of energy exactly determines the scal- 
ing of the third order structure function of velocity field 
which is responsible for carrying the energy flux. This 
statement is known as the Kolmogorov |-law. It states 
that for length scales, r, in the inertial range, the third 
order structure function is: 

{{ Vl (r)-v l (0)] 3 )) = ~er, (1) 

where e is the energy flux, and m is the longitudinal com- 
ponent of the velocity. The above law does not depend on 
the dimensionality of space or details of energy injection 
and dissipation. This is an important benchmark in the 
theory of turbulence for which no controlled approxima- 
tion schemes are known for the computation of general 
correlation functions. Analogous results have been ob- 
tained for other hydrodynamic turbulent systems, such 
as magneto-hydrodynamics [3, Burgers turbulence [li| 
and turbulent advection (see [20j and references within) . 
If, based on the |-law, one postulates, as Kolmogorov did 
in his 1941 theory, that all velocity correlation functions 
depend only on the energy flux as a parameter, then one 
is lead to predict that the n th order structure function 
should scale with exponent ^. In fact, experiments in- 
dicate that turbulence exhibits multi-scaling: the scaling 
exponent of the n th order structure function varies non- 
lincarly with n. As was realised later by Landau, Kol- 
mogorov's 1941 theory is an example of 'mean field' the- 
ory which is valid only if energy flux fluctuations around 
the mean value are small. 

The relevance of this discussion of universal features of 
turbulence for the present paper, ostensibly on a rather 
different topic, is outlined in a recent series of papers 
[Tol [TTI . I2H ] in which we explore the implications of the 
analogy between energy cascades in turbulence and mass 
cascades in the Takayasu model. We make extensive use, 
in the present work, of the universal scaling properties of 
the flux-carrying correlation function as expressed in the 
CFR [llj - the direct analogue for the Takayasu model 
of the |-Law. For a more detailed discussion of the CFR 
in the context of cluster aggregation see [22j . 

The rest of the paper is organised as follows. We first 
(Sec. |n|) define the model and briefly summarise known 



results on the phase transition and the behaviour of cor- 
relation functions in different phases. We then (Sec. IIII|) 
use the first equation in the BBGKY hierarchy (Hopf 
equation) relating the one- and two-point functions of 
the model to prove the existence of the growing phase of 
the model and to establish a rigorous upper bound on 
the position of the transition point. The stationary Hopf 
equation (at t — 00) is then re-expressed fSec.HVp as the 
balance condition between processes of aggregation, de- 
position and evaporation. The implications of this bal- 
ance condition are then analysed in the growing phase 
(Sec. E} and in the exponential phase (Sec. IVI|) . The 
existence of an asymptotically constant flux of mass in 
the growing phase is established and used, in conjunction 
with the CFR, to calculate the universal scaling of the 
two-point function in all spatial dimensions. In the expo- 
nential phase, the flux decays exponentially and, using a 
combination of mean field approximation and numerical 
methods, we show that the CFR is replaced by a scale- 
by-scale balance between mass fluxes due to evaporation 
and aggregation. Finally, (Sec. IVII|) the scaling theory of 
the model near the critical point is developed, again mak- 
ing use of the CFR to obtain a relation between scaling 
exponents, and applied to the calculation of correlation 
functions at the critical point. 



II. THE MODEL 

Consider a d-dimcnsional hyper-cubic lattice. A lattice 
site may be occupied by at most one particle. Let the 
mass of the particle at site i be denoted by rrii, where 
m,i = 0,1,2,.... By convention, rrii = corresponds to 
an empty site. Starting from an initial state where all 
the lattice sites are empty, the system evolves in time 
through the following Markovian stochastic processes: 

1. Diffusion-coagulation: With rate D, a particle hops 
to a nearest neighbour site. If the target site is 
already occupied, then the two particles coagulate 
to form one particle whose mass is the sum of the 
two particles. If i is the site and j the target site, 
then mi — > 0; nij —> mj + m^. 

2. Deposition: With rate q, a particle of mass 1 is 
adsorbed at a site i, increasing the mass at the site 
by one, i.e., m t — > rrij + 1. 

3. Evaporation: With rate p the mass at a site is de- 
creased by one, provided there is a nonzero mass at 
that site, i.e., mi — > mi — 1 + 5 mi fl. 

By an appropriate renormalization of the time scale, the 
diffusion rate D is set equal to 1. The model is then 
characterised by two independent parameters: p - the 
evaporation rate and q - the deposition rate. We refer to 
the model as the evaporation-deposition model or EDM. 

Let P(m, t) be the probability that a site is occupied 
by a particle of mass m at time t. More generally, let 
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P(rrii 1 , rrii 2 , . . . , mj„ ; t) be the the n-point joint probabil- 
ity distribution that lattice sites i±, 12, • ■ • , i n are occupied 
by particles with masses , rrii 2 , . . . at time t. In 
the limit of large time, the EDM approaches a steady 
state. All probability distribution functions accordingly 
become independent of time. We will refer to the proba- 
bility distribution functions in the steady state by drop- 
ping the time argument, e. g., P(m) = lim^oo P(m,t). 

A striking feature of the EDM is the existence of a non- 
equilibrium phase transition at t = 00 in all dimensions. 
This was established in [l3j using a combination of mean 
field analysis and numerical simulations. For fixed evap- 
oration rate p and sufficiently small deposition rate q the 
steady state has finite total mass (m) = J2m=i mP(m). 
If however, q is increased beyond a certain critical value, 
q c (p), the total mass in the steady state becomes infinite. 
It has been also observed, that for q < q c (p), the large 
mass asymptotics of P(jn) are exponential, whereas for 
Q = Qc(p) and q > q c (j>) the asymptotics of the mass 
distributions are algebraic but with different exponents: 




raj 'm* 



when q < q c , 
when q = q c , 
when q > q c , 



(2) 



where m* is a q dependent cut-off, and r 7^ r c , see [13j, 
IMEl for more details. The three phases of the EDM at 
t = 00 will be called the exponential phase (q < q c ), the 
critical phase (q = q c ) and the growing phase (q > q c ). 



III. THE PHASE TRANSITION 

EDM is a Markov model with configuration space Z^ d . 
Accordingly, the evolution of the probability measure on 
the configuration space is governed by the master equa- 
tion. Marginalizing the master equation over all but one 
lattice site, it is straightforward to derive the system of 
equation describing the time evolution for P(m,t): 



dP(m, t) 
dt 



= q[P(m-l,t)-P(m,t)] (3) 
+ p[P(m+l,t) - P{m,t)\ 

OO 

- V] P(m,m',t) - P(m,t) 



m' — 1 
rn — 1 



dP{0,t) 
dt 



P(m f , m — m! , t), m > 0, 

re'=0 

-qP(0)+pP(l) (4) 

OO OO 

P(0,m') + ^2 P( 



E 

ra' — 1 



where P(mi,rri2,t) is the joint probability distribution 
function of two neighbouring sites having mass mi and 
Tri'2 . The simplest way to verify and ([5]) is to enu- 
merate the number of ways a site can be occupied by a 



particle of mass m and the number of ways a site already 
occupied by a particle of mass m can change to another 
mass in a small time dt. These equations arc not closed: 
it is the evolution equation for the one-point distribution 
of mass, but its right hand side contains the two-point 
joint distribution function. In fact, (j4]) and ((5]) are the 
first members of the infinite BBGKY (Hopf) hierarchy, 
which in general is impossible to solve without resorting 
to closure approximations such as mean field theory or 
perturbative renormalization group. 

There is, nevertheless, a great deal of information 
about phases of EDM which can be extracted from the 
first members of the hierarchy without resorting to any 
closure assumptions. For example, in the growing phase 
of the EDM, there is a hidden conserved quantity - the 
flux of mass due to aggregation, which will allow us to ex- 
tract the exact large mass asymptotics of the two-point 
mass distribution from ^) and ([5]) in all space dimen- 
sions. 

We start by demonstrating the existence of the grow- 
ing phase of EDM in all dimensions without resorting 
to any closure approximations such as mean field theory. 
Multiplying both sides of (U) with m and integrating over 
all masses, it is easy to find that J(t), the flux of mass 
into the system, is given by 



TM _ d(m(t)) 
J(t) = — - — =q-ps(t), 



(5) 



where s(t) = Ylm=i P{ m i t) is the probability that a site 
is not empty and (m(t)} is the average mass at time t. 
As s(t) < 1, 



J{t)>q-p. 



(6) 



Therefore, if the deposition rate q is larger than the evap- 
oration rate p, then the flux of mass into the system is 
positive at all times, implying that the system is in the 
growing phase with (m(oo)) = 00. Equation (J6]) provides 
a simple upper bound for the critical value: q c (p) < P- 

The above upper bound can be improved quite easily. 
Equation ([5]) for m = may be rewritten as 



dP(0,t) 
dt 



1 - (2 + q)P(0, t) +pP(l, t) + P(0, 0, t). (7) 



Dropping the non-negative terms pP(l,t) and P(0, 0,i), 
we obtain the inequality 



dP(0,t) 
Jt 



>l-(2 + «)P(0,t), 



(8) 



for all times t > 0, with the initial condition P(0,t = 
0) = 1. This reduces to the following bound on P(0,t): 
P(0,t) > F(t), where F(t) is the solution of the linear 
differential equation 



dF(t) 
dt 



= l-(q + 2)F(t), P(0) = 1. 



Solving the differential equation, we conclude that 



P(0,t) > 



1 



f 



"(9+2)t 



(9) 



(10) 
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The derived inequality immediately leads to an upper 
bound on particle density s = 1 — P(0, t = oo): 



In the exponential phase the global balance condition ([5]) 
dictates that s = q/p which is consistent with the bound 
(HH only if 

q < <H B = \ (p ~ 2 + Vp 2 +4) , (12) 

which establishes an upper bound on q c (p) that is tighter 
than the simple q c (p) < P bound. 

Proving the existence of the exponential phase or 
equivalently deriving a lower bound for q c (p) is difficult. 
However, numerical simulations and a non-trivial mean 
field analysis carried out in [l3[ shows that as we decrease 
evaporation rate q while keeping p constant, the system 
undergoes a non-equilibrium phase transition into an ex- 
ponential phase in which the total mass present in the 
system remains finite at all times. The transition occurs 
at some critical value of the deposition rate q = q c {p) > 0. 
Within the mean field approximation [l3j 

q™ F (p)=p + 2-2y/p-+l. (13) 

Numerics [see figured] suggest that the curve q = q^ IF {p) 
is a lower bound on the true phase transition bound- 
ary. We cannot prove this conjecture and thus establish 
the existence of the exponential phase by using just the 
first pair of equations from the BBGKY (Hopf) hierar- 
chy. The relationship between the rigorous upper bound 
(fl"2|) , the mean field estimate of the phase boundary (fT3|) 
and the numerical estimate of q c {p) in one dimension is 
illustrated in figure [TJ We now discuss the large time 
statistics of the model in each of its phases. 

IV. MASS BALANCE IN THE STEADY STATE. 




0.5 1 1.5 2 2.5 3 
Evaporation rate, p 



FIG. 1: Phases of evaporation-deposition model. Squares: 
numerical measurement of q c (p) in one dimension; Dashed 
curve: mean field result ()13|) ; Solid curve: upper bound (|12l) : 
circles: Monte-Carlo simulation of one-dimensional EDM 
with long range hopping. 
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FIG. 2: Balancing fluxes in the steady state of the model. 



In the steady state the amount of mass contained in 
every finite interval [l,m] is constant: 



dt 



5>P(/i) = 0. 



(14) 



Using (j4]) this can be re- written as the following balance 
condition: 



where 



j i j( m ) — r(i) i j( m ) 

J-m ~ < J aqq ° aqq ' ° ev ) 



(i=i fi=i 



(15) 



(16) 



is the bulk term equal to the rate at which mass escapes 
from the system due to evaporation and deposition, 



411 = <? p (°)> 



(17) 



is the flux of mass into the interval [l,m] due to deposi- 
tion of particles of unit mass onto empty sites, 



jM = pmP(m + 1) 



(18) 



is the flux of mass into the interval [1, m] due to evapo- 
ration of particles of mass m + 1, and finally 



fj,=0 fj,'=0 



(19) 

is the flux of mass out of the interval [l,m] due to 
diffusion-aggregation. Figure [2] illustrates the balance 
equation graphically. The balance (|15[) is realized differ- 
ently in different phases of the model, which is elaborated 
upon below. 
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FIG. 3: Numerical measurement of aggregation flux in the 
growing and critical phases of EDM in one dimension. The 
simulations are for p = 1.0, when q c m 0.3072. For q > 
q c (p), J a gg(m) m ^°° J > 0. At the critical point, the large 
mass limit of J a „„ is zero. 



the Zakharov transform [IT| to obtain 
P(mi,m 2 ) = 7 rrj^f ( — ) , mi,m 2 > 1, (24) 

where f(x) = f{l/x) is an unknown scaling function. 
This result holds in all dimensions and is derived with- 
out any closure assumptions. We conclude that the two 
point function in the growing phase of the EDM scales 
exactly as in the aggregation model with deposition and 
no evaporation. 

It is worth noting an alternative representation of the 
above result in Fourier space. Let 

oo 

iKz!,Z2)= £ P(m u m 2 )(z^ -l)(zp -1). (25) 

TJl 1,2—0 

Then the equivalent form of f[24|) in terms of ip(zx, Z2) is 
^z,z) = -J(z-\) + 0[{z-\)\ (26) 



V. GROWING PHASE 

Consider the large mass limit of (|15|) . Being a proba- 
bility distribution, P{m) sums to one, which means that 
P(m) must go to zero faster than 1/m as mass m tends 
to infinity. Therefore, ji™^ vanishes for large masses, 

jM = pmP( m + 1) '™ 0. (20) 

The bulk evaporation term I m approaches a limit: 

I m m ^?(p-q)s, (21) 

where s = 1 — P(0). Therefore the balance equation 
acquires the form 

4$ = 4H -Ioo + 0(m- a ) = J + 0(m- a ), Vm » 1, 

(22) 

where a > and J is the total flux of mass at m = oo, 

J = q-ps, (23) 

as in ([S]). In the growing phase J > and we find that 
the flux of mass due to diffusion-aggregation is asymp- 
totically constant in the limit of large mass despite evap- 
oration and deposition being present at all scales. In 
figure [3l the variation of Jagg with m is shown for differ- 
ent values of q > q c . For large m, the current goes to a 
non-zero constant in the growing phase. 

We conclude that the large mass limit of J^J in EDM 
is identical as that of the model with no evaporation, but 
with a renormalised deposition rate J = q — ps. The fact 
that there is a constant flux of mass in the model allows 
us to apply the constant flux relation derived in [ll| to 
determine the exact scaling of the two-point function in 
the growing phase of the model. In particular, the equa- 
tion JagJ = constant for large m may be solved using 



which is strongly reminiscent of the Kolmogorov's 4/5- 
th law in fluid turbulence. The easiest way to derive 
(f2"6"| is to differentiate the generating functional ip(z) = 
^^dPfmjz" 1 with respect to time and resolve the re- 
sulting time derivative using the equations of motion Q . 

Equation (|24|) is a direct consequence of the asymptotic 
constant flux condition ([2^]) which holds in the growing 
phase despite the presence of evaporation. It is natural to 
conjecture that the growing phase of EDM is equivalent 
to aggregation model with suitably renormalised deposi- 
tion rate. A consequence of such a conjecture would be 
the expression for all scaling exponents for multi-point 
probability distribution functions of the model in one di- 
mension. Namely, based on the results of (lol | , we expect 
that 

P(Atoi, Am 2 , . . . , Am n ) — A~ 7 "P(mi, m 2 , . . . , m n ), 

(27) 

where in one dimension 

An n(n — 1) , . . 

In = y + V 6 , rc= 1,2,3,... d=l (28) 

We check this conjecture numerically. Let 

P fe (m) = P{ m,m,...,m) , (29) 

k 

where the masses m are on neighbouring sites. We mea- 
sure Pfc(m) in the growing phase in one dimension. The 
results arc shown in figure |4l The data are in good agree- 
ment with (|28jl . A rigorous explanation of the observed 
result that the growing phase of the aggregation model 
with evaporation rate p > and deposition rate q > is 
equivalent to the aggregation model with no evaporation 
and deposition rate q' = J = q — ps remains an open 
problem. 
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ln(m) 

FIG. 4: Numerical verification of ([27)l and ([280 for one-, two- 
and three-point correlation functions in one dimension. The 
theoretical values of the scaling exponents are given by (|28|) : 
7i = |, 72 = 3, 73 = 5. 
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FIG. 5: Numerical measurement of currents in the exponen- 
tial phase for one-dimensional EDM. The aggregation cur- 
rent J ag g(m) and mass distribution P(m) decay exponentially 
with mass in quantitative agreement with mean field theory. 
Inset: Numerical verification of (|36[) in one dimension. 



VI. EXPONENTIAL PHASE 

We now analyse the implications of the balance equa- 
tion (|T5|) in the exponential phase. In this case the total 
flux of mass at infinite mass is zero. As a consequence, 
we claim that the large-mass behaviour of the model in 
the exponential balance is described by a steady state 
characterised by a scale-by-scale balance between aggre- 
gation and evaporation fluxes. As far as we are aware, the 
non-perturbative analysis of the exponential has never 
been attempted. Therefore, we will base our arguments 
mostly on mean field results and numerical simulations 
in one dimension. 

The mass fluxes , I m and Jagg defined in (fT5|) , (TP?)) 
and (TP5j) can be written for large masses as 



Ji m) ra pmP(m) 



(30) 



and 



and 



J (1) 

199 



A(p-q)- 



7.3/2 



-m/m* 



(34) 



Thus, Jagg—Im decays with mass faster than the flux due 
to evaporation. The balance equation (|I5p then reduces 
to 



j{m) r(m) ^ i 
°agg ~ " ev i ^ L - 



(35) 



J&l-Itn « (q-ps)-(q-p) / dflP(fl) = (p-q) / d/iP(/i) 

J rn J m 

(31) 

where the global mass flux q — ps has been set to zero 
for the exponential phase. Clearly, if P{m) goes to zero 
at large masses faster than a power law, for example as a 
stretched exponential or as an exponential, then ji™' 1 S> 
J (1) - I 

u agg * un - 
it has been shown in that, within the mean field 
theory, for large masses, P(m) has the form 



Therefore, a scalc-by-scale balance is established in the 
exponential phase for each value of m. The flux of mass 
out of the interval [l,m] due to diffusion- aggregation is 
equal to the flux of mass into the interval [1 , m] due to 
evaporation of particles of mass m + 1. Simultaneously 
with this scale-by-scale balance, the global balance holds: 
the rate of the total mass increase due to deposition is 
exactly equal to the total rate of mass loss due to evap- 
oration. Using the explicit expression (|T8|) for it 
is possible to re-write ([33]) as a relation between mass 
distribution and the aggregation flux: 



P(m+ 1) 



1 

pm 



r(m) 
agg I 



(36) 



P(m) 



A 



7,3/2 



e"™/™*, m»l, 



(32) 



where m* is the effective mass determined by the singu- 
larities of the mean field moment generating function. It 
follows that 



J 



pm~ 1 l 2 e- m ' m ' 



(33) 



which replaces the CFR of the growing phase. 

Our conclusion that in the absence of global flux, the 
steady state is described by a scale-by-scale balance be- 
tween evaporation and aggregation seems natural and 
should hold irrespective of the validity of the mean field 
assumption. The exponential nature of P(m) should be 
true in one dimension too, though corrections to it would 
be different from the mean field result. In figure [5j we 
show the results of simulation in one dimension. The 
data are consistent with (f3"6")l . 
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VII. SCALING THEORY OF THE CRITICAL 
PHASE 

In this section, we derive the scaling exponent for the 
two-point probability distribution function P(mi, 7712) at 
the critical point q = q c (p) in terms of the scaling ex- 
ponents of the one-point distribution P(m). Theoretical 
analysis and numerical simulations show that the aggre- 
gation flux vanishes as m — > 00 at the critical point, see 
figure [3] Therefore, the calculation of the two-point dis- 
tribution based on constant flux carried out in Section 
|V]does not apply to the critical point. Instead, we will 
base our argument on the exact result (f2~4"f valid just 
above the critical point and the scaling theory for the 
EDM developed in [HHI. 

Let us fix the evaporation rate p. The parameter Sq = 
q — q c (p) > measures the distance to the critical point 
in the growing phase. We are interested in the limit of 
Sq — > +0 and t — > 00. In this limit P(m, q, t) displays the 
scaling form [l3|, [23[ 



P(m, Sq, t) 



n T <= V 



in 



(37) 



in terms of three unknown exponents 0, a, t c , and the 
scaling function of two variables Y\. Of interest is one 
more exponent 9 defined as follows. In the growing phase, 
the average mass increases linearly with time, and the 
mean rate of increase scales with distance to the critical 
point as J ~ Sq 6 , where 9 > is the growth exponent. 
In other words, 



(m) ~ Sq e t, 9>0. 



(38) 



The growth exponent 9 is not independent of the expo- 
nents r c , <p and a which define the scaling of mass dis- 
tribution: calculating the expected mass (m) using 
and comparing the answer with (|38[) we find that [L; 



' [1 - a(2 - t c )] 



(39) 



The exponents a, </> and r c are not independent, and 
are connected through the scaling relation, [24j : 



a(r c 



D = *. 
' 2 



(40) 



The derivation of (|4T)]) is not as straightforward as (|3"9")) 
and is based on a catchment area argument, see [24| for 
details. We conclude that in the scaling regime, the mass 
distribution is characterised by two independent scaling 
exponents, e.g. a and <p. 

Now let us consider the two point probability distri- 
bution. In the limit t — > 00, it will have the following 
scaling form 



P(m, m, Sq) ~ — jt^2 (mSq^) 

TKl tic ^ ' 



(41) 



where Y2 is a scaling function of one argument. Note 
that (|41[) is based on the assumption that exponent <p 



determines the behaviour of all probability distribution 
functions near the critical point. In the growing phase, 
the scaling of P(m, m, t) as a function of mass m and 
flux J is determined by the constant flux relation: 



P(m, m, t) 



J 



m > 1, 



(42) 



where h = 3. 

If Y 2 (x) ~. 
(|42|) leads to a pair of equations 



M for x 3> 1, the comparison of (|4T|) and 
" for exponents (i and h c : 



fx = h — h c 



Solving for h c we find that 



1 



h c = - [d - 2 + r c (2 + d)} , 



(43) 



(44) 



which is consistent with conclusions of [14| . We conclude 
that at the critical point 

/ x x (2-d)-T c (2+d) 

P(m,m) | g =g (p)~ [m^J . (45) 

Note that at the critical point the scaling of the two- 
point function does depend on d, the dimensionality of 
space. This is to be contrasted with the universal scal- 
ing law of P(m,m,t = 00) in the growing phase due to 
the constant positive flux of mass in mass space. Note 
also the independence of h c on the evaporation rate p. 
In figure [51 we measure P(m, m) in one dimension at the 
critical point. The exponent r c in (|45|) is unknown, but 
can be eliminated in terms of the mass distribution P(m). 
The numerical results in figure [6] confirm the scaling law 
(1451) . It is quite surprising and encouraging that the con- 
stant flux relation valid in the growing phase combined 
with a simple heuristic scaling theory allows one to make 
correct guesses concerning the statistics of the model at 
the critical point. A rigorous derivation of (|45[) remains 
an open problem. 



VIII. SUMMARY AND CONCLUSIONS 

To summarise, we have studied the phase transition 
between a growing phase (in which the total mass grows 
for all time) and an exponential phase (in which the total 
mass is bounded) in a lattice based diffusion-aggregation 
model with injection and evaporation. This transition 
is controlled by the ratio of the rate of deposition of 
monomers, q and the rate of evaporation, p. We es- 
tablished a rigorous upper bound on the transition point 
q c {p) valid in all dimensions and showed that the growing 
phase has an asymptotically constant flux of mass with 
an associated CFR for the two-point probability distri- 
bution function, again valid in all dimensions. From a 
combination of theoretical analysis and numerical simu- 
lations, we conjectured that the growing phase is in the 
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FIG. 6: Numerical verification of (|45[) in one dimension: For 
any r c , g5| and j37J imply that for d = 1, ~ m°. 

Inset: P(m,m) scales as m -4,44 . This is consistent with /i c = 
| corresponding to the value of r c = ^ argued for in [23 ], 



universality class of the aggregation model with deposi- 
tion but no evaporation, the only effect of evaporation 
being the modification of the value of the mass flux. Nu- 
merical simulations in one dimension show that all prob- 
ability distribution functions of the model in the growing 
phase scale with mass exactly as they do in the model 
without evaporation. It would be interesting to verify 
this theoretically. 

The existence of a non-trivial exponential phase, (i.e. 
the fact that q c (p) > 0) was only established at the 
level of mean field theory although numerical evidence 
for the existence of this phase in one dimension is in- 
controvertible. The presence of the phase transition can 
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